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1 General equation for the calculation of amplitudes 

It is well known that the high order calculations of the observables within the perturba- 
tion theory has the serious difficulties (especially, if we take into account the polarization 
effects). Method of the direct calculation of amplitudes simplifies these calculations es- 
sentially. 

There is an even number (2A^) of fermions in initial and final state for any reaction 
with Dirac particles. Therefore every diagram contains N non-closed fermion lines. The 
expression 

Mif = UfQui = TT{QuiUf) (1.1) 

corresponds to every line in the amplitude of the process, Ui, Uf are the Dirac bispinors 
for free particles and Q is a matrix operator characterizing the interaction. The operator 
Q is expressed as a linear combination of the products of the Dirac 7-matrices (or their 
contractions with four-vectors) and can have any number of free Lorentz indexes. 
We must obtain the expression for the operator 

UiUf (1.2) 

in JTH). 

Note that Dirac equation 

pu{p,n) = mu{p,n) , (1.3) 
and equation for the axis of the spin projections 

-~^^hu{j),n) = ±u{p,n) (1-4) 
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as well as normalization condition for bispinor define u{p,n) up to the phase factor {p 
is four-momentum, n is the four-vector specifying the axis of the spin projections and 
o = IfiCL^ for any four- vector a) ^. 

Let use this fact to construct the operator ()1.2|) : 

UiC::^Ui— — I = . =u , (1.5) 



\UiU\ 



Tr(PP,) 



where u is an arbitrary bispinor and V are projection operators for Dirac particles. 

V{p, n) = u{p, n)u{p, n) = -{p + m)(l + 7^n) (1.7) 
is the projection operator for a particle with mass m, 

p^ = , = —1 , pn = , UU = 2m . 
For massless particle the projection operator has form: 

u±{p)u±{p) = ^(1 ± l')P = 'P±ip) (1-8) 

where: 

p'^ = , u±'y^u± = 2p^ 

(signs ± correspond to the helicity of particle) . 
As a result we have 

UiUf ~ , . (1.9) 

^TiiVPi) Tr{VPf) 

It follows from (jl.Sp . ()1.6|) that both Ui and uj obtain individual phase factors, i.e. direct 
as well as exchange diagrams multiplied by the same phase factor, and consequently one 
can ignore it. Therefore we will use the equality sign instead of the symbol ~ in formulae 
for amplitudes. 

Finally, we have the general formula for the calculation of amplitudes 

,,Q„.=^(22m^. (1.10) 

/Tr(TO) TiiVVf) 



^We use the same metric as in the book Q]: 

= (ao, a), = (ao, -a), ab = a^6^ = ao^o - ab, 7^ = i7°7^7^7^ 
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This expression enables to calculate the amplitude numerically. The complex values 
obtained in this way may be used to evaluate the squared amplitude. 

Note that this method can be generalized easily to the case of the reaction with 
antiparticles. To do this it is sufficient in ()1.1U|) to replace the projection operators of 
particles by the operators of antiparticles. 

More details may be found in j2]. 



2 Example: the calculation of electron-electron scat- 
tering amplitudes 

As an illustration of application of the method we will consider the lowest order amplitudes 
of electron-electron scattering. Two Feynman diagrams presented in Fig^ correspond to 
this process. 



Pi ' r > P^ Pi ' r ' ■^'4 

rPi - P3 \Pi - Pa 

P2 ' J > Pa P2 > J ' P3 

1 2 

Figure 1: Diagrams for electron-electron scattering 

^ (^p -p )2 ^3^P^1 ■ ^47''^2 , (2.1) 
Ms = --. ^ UAlpUl ■ U3Y'^2 ■ (2.2) 

There are two variants of formula (jl.lUj) for massive Dirac particles: 

_ Tt[{1 - -f^)q{pf + mhf + m - pfhf)Q{pi + mfii + m + pihi)] 

8\J{qPi) + m{qni)^{qpf) + m{qnf) 

if P = -(1 — 7^)g , (the massless four- vector q can be chosen arbitrary, but it has to be 
the same for all non-closed fermion lines of the considered diagrams), and 

_ Tr [(1 + r')q{pf -mnf + m + Pfnf)Q{pi - mfij + m - pifij)] 
{(IPi) - m{qni)^{qpf) - m{qnf) 
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if p = -(1 + 7^)^ . Relation between (Q and gives by 



UfQui = (ufQui)' 



- Tr [(1 + 'y^)qpihiqpfhf] 



(2.5) 



For definiteness we will use fl2.H|l . 

Note, to check formula ()2.5|1 it necessary to use identity: 

(1 ± 7')gQ(l ± 7^)^ = Tr [(1 ± ^')qQ] (1 ± ^')q , 



(2.6) 



where an operator Q is any matrix. (The proof of ()2.6|) is given in pi.) Besides, we will 
use well-known formulae of 7-matrix algebra: 



7^g2n+i^p = -2Q 



R 



and 



1 ± 7')7pQ'V = (1 ± 7') Tr (1 T l')Q 



(2.7) 



(2.8) 



where Q^""*"^ is a product of odd number of the 7-matrices; Q/f"^^ is a product of the 
same 7-matrices, rewritten in reversal order and is a product of even number of the 
7-matrices. 

Thus, we obtain from ()2.3p : 

_ Tr [(1 - 7^)g(p3 + mfi3 + m -^3^3)7^(^1 + mfii + m + pifii)] 
UsIpUi — — ^ . , n . ' l^-^/ 



M47''W2 



Tr [(1 — 7^)g(p4 + mn^ + m — P4n4)'-f^{p2 + mn2 + m + ^2^2)] 



{qP2) + rn{qn2)^J (qpi) + m{qni) 

1 



(2.10) 



V (^^'i) + '^iQ^i)\JilP2) + m{qn2)^J{qp3) + m^qn^i)^ {qp^) + m(gn4 



|Tr (1 - 7^)g(j03 + m?T,3)7p(j3i + mni) Tr (1 - 7^)g(j94 + mn4)7^(j32 + "^^2) 



+ Tr 



(1 - 'y^)q{p3 + mh3)^p{pi + mni) Tr (1 - Y')q{m - p4h4)Y{m + ^2^2) 



+ Tr (1 - 7^)g(m - P3^3)7p('"^ + Pi^i) Tr (1 - 7^)g(p4 + mh4)Y{p2 + "^^2)] 

+ Tr (1 - 7^)^(m -^3^3)7^(777, + pini) Tr (1-7^)^(771-^477,4)7^(777 + ^2^2) } 

(2. 
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Using (j2.6j) . (j2.7j) . for first term in (j'i.llj) we have: 



Tr 



(1 - 7^)g(p3 + "^^3)7p(Pi + rnhi) Tr (1 - 7^)g(p4 + mhi)Y{P2 + ^^2) 



2Tr 



(1 - 7^)g(p3 + "^^3)7p(Pi + mhi)q{p4 + mhi)Y{P2 + mn2) (2.12) 



-4Tr 



(1 - 7^)^(p3 + 'mhs){p4: + ■mfii)q{pi + mhi){p2 + mfi2) 



For second term in (j2.1ip we obtain through (j2.6j) . (j2.8j) 



Tr 



(1 - 7^)g(p3 + mn3)7p(pi + mni) Tr (1 - 7^)g(m - p4n4)7''(m + p2'"'2) 



2Tr 



(1 - 7^)g(p3 + mns)-fp{pi + mhi)q{m - pihi)Y{m + ^2^2) (2.13) 



2Tr 



(1 - 7^)^(p3 + mn3)(m + ^2^2) Tr {I - ^^)q{m - pihi){pi + mfii) 



In the same way 



Tr 



(1 - 7^)g(m -p3n3)7p(m + ^1^1) Tr (1 - 7^)5(^4 + mn4)7''(p2 + "^^^2) 



2Tr 



{I - '^^)q{pi + mhi){m + pihi) Tr (1 - 7^)^(m - p3n3)(p2 + "^"^2) 



and 



Tr 



(1 - 7^)g(m - P3n3)7p(m + pini) Tr (1 - 7^)g(m - P4n4)7''(m + ^2^2) 



= -4Tr 
As a result we have 



(1 - 7^)g(m - p3ns){m + p4hi)q{m - pihi){m + P2'^2) 



(2.14) 



(2.15) 



32^^ (qpi) + m{qni)^{qp2) + m{qn2) {qpz) + m{qn^)^{qpi) + m(gn4) 
I - 2 Tr (1 - 7^)g(j53 + mhs){p4 + mn4)q{pi + mhi){p2 + mn2) 
+ Tr (1 - 7^)g(p3 + mr2.3)(m + ^2^2) Tr (1 - 7^)g(m - p4n4)(pi + mni) 



+ Tr 



(1 - 7^)g(p4 + mn4)(m + pifii) Tr (1 - 7^)g(m - P3fi3)(p2 + "^^^2) 



-2Tr 



(1 -7^)g(m-p3n3)(m + p4n4)g(m-pini)(m + j32n2) | 

(2.16) 



In principle we may use (j2.16|) for calculation of amplitude but it possible to simplify 
essentially the last term in this expression. Really, using consequence of formula (j2.5|) 



Tr 



Tr 



(1 - 7^)g(p2 - mn2)Y{p4: - mhi) 



Tr [(1 - 7^)gP4^4#2?^2] 



8 [{qP2) - m{qn2)] ■ [{qpi) - m{qn4)] 



we have instead of ()2.15|1 



(2.17) 



Tr 



(1 -7^)^(m-p3n3)7p(m + pini) Tr {1 - j^)q{m - p4h4)Y{m + P2n2) 



-2Tr 



(1 — 7^)g(p2 — 'T^?^2)(^ + Pini) Tr (1 — 7^)^(m — p3n3)(p4 — mn4) 

Tr [(1 - 7^)gP4^4#2^2] 



(2.18) 



8 [{qp2) - m{qn2)] ■ [{qpA) - miqu^)] 
Similarly, for first term of ()2.16|) . using the other consequence of formula ()2.5p 



Tr 



1^1 - 7^)g(p4 + mn4)Y{p2 + mn2] 



Tr 



1 - Y)q{'m + p2h2)Y{m - ^4^4) 



Tr [(1 - 7^)gp4^4#2^2] 



8 [iqP2) - m{qn2)] ■ [{qpi) - miqn^)] 



we obtain instead of (j2.12|) 



(2.19) 



Tr 



(1 - 7^)g(p3 + mn3)7p(^i + mni) Tr (1 - 7^)g(p4 + 772^4)7^(^2 + "^^2) 



-2Tr 



(1 - 7^)g(p3 + mn3)(?77 - ^4774) Tr (1 - 7^)^(777 + p2^2)(Pi + "^fii) 

Tr [(1 - Y)QPi'f^iQP2n2\ 



(2.20) 



8 [{,qP2) - m{qn2)] ■ [{qpi) - m{qni) 
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The final equation is 



UalpUi ■ U4YU2 



■|Tr (1 - 7^)g(p3 + mn3)(m + ^2^2) Tr (1 - 7^)g(m - j94?7,4)(pi + mui) 
+ Tr (1 — 7^)^(p4 + 771,^4) (m + pifii) Tr {\ — '^^)q{m — p^n^){j)2 + 'fTin2) 
— {Tr (1 — 7^)^(p3 + mn3)(m — J94n4) Tr {1 — '~^^)q{m + p2n2){pi + mfii) 
+ Tr {1 — '^^)q{p2 — mn2){m + pifii) Tr (1 — 7^)^(m — P3n3)(p4 — 771^4) } 

Tr [(1 - 7^)gp4"'4#2^2] 



}■ 



8 [{,qP2) - m{qn2)] ■ [{qpi) - 777(5774)] 

(2.21) 

The expression for 

is obtained from ()2.2H) by replacing the subscripts 3 4. 

The obtained formulae are simple enough and may be used for numerical calculation 
of the amplitudes. 

If the denominators in ()2.2H1 become zero for some values of the vectors of problem, 
then one needs to change the values of arbitrary four- vector q in all formulae. 



3 Construction of polarization vectors of photons from 
the vectors of the problem 

If photons participate in reaction, we must to construct their polarization vectors from 
other vectors coming in amplitude calculation. It is necessary to use Gram - Schmidt 
orthonormalization process to do this (see j^). 

As it shown in |3], in the Minkowski space, this process give 

1. Let p is an arbitrary four- momentum, such that 

= 777^ 7^ , 

a, h and c are arbitrary vectors. Then four vectors Iq, li, I2, I3 form an orthonormal 
basis: 

/o = -, (3.1) 



7 



G 



p a 



m^ap - {pa)pp 



m 



-G 



p a 



.P OL 

Gi^ a b 

.pap 



.1/2 



mJ {pay — rri^a^ 



(3.2) 



pa) \p a b 



1/2 



[a (pb) — {pa){ab)] Pp + [m [ab) — {pa){pb)] ap + [{pa) — m a ] bp 
^ {pay — m?a'^\j2{jpa){pb){ab) + m?a'^b'^ — m?{abY — a^^pb)"^ — h^^paY 



G 



{I 



3jp 



p a b c 
p a b p 



(3.3) 



Q I P a ^\ Q f P ^ ^ ^ 

p a bj \p a b c 



1/2 



where G are Gram determinants. Finally we take into account that in the Minkowski 
space 

' p a b 



G 



p a b c 



p'^a^'Wee^px^p'^a^b^'' 



Gil \ '')=ep.„rP^a'b-c^ep^^xp''a^b\ 
.P a p J 



and last vector becomes: 

^pa(i\P'^ 



raf^b^ 



3)p 



£pal3Xp 



G 



p a b 
p a b 



^2{pa){pb){ab) + m'^a^b'^ — m?{abY — a'^{pbY — b'^{paY 

(3.4) 



Note that in the final formulae the vector c is absent. 

So, using the three vectors (one of which is a four-momentum of a massive particle) 
and a total antisymmetric Levi-Civita tensor, an orthonormal basis can be always 
constructed in the Minkowski space. 

2. Let us consider now the construction of the basis for reaction with massless particles. 
Let ki and k2 are four-momentums of the problem, such that 

• 

We consider a particular form of the basis ()3.1|) - ()3.4|) at 

p = ki + k2, m= \j2{kik2) , a = k2 : 



(^2)p = 



ki + k2 
h — I I 
V2(A:ifc2) 

^ _ -ki + A:2 

V2(A;iA;2) ' 

-(fc2b)(A;i),-(fcib)(A;2),+ (A;ifc2)bp 
y^2(A;iA;2)(A;i6)(fc26)-fo2(A;iA;2)2 

^ pcxi3\k'ik2h^ 

Sj2{k,k2){k,b){k2b)-b^{kik2f ' 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



If ki is four-momentum of photon, polarizations vectors for it may be constructed from 
()3.7j) . ()3.8|1 . where it is appropriate to choose b = q, [q is the same massless vector which 
contained in expressions for amphtudes ()2.3p . ()2.4|) ]. 



-(fc2g)fcip + {kik2)qp - {kiq)k2p 
2{kik2)ikiq){k2q) 



(1 



3)p 



^papxkiq^k. 



2{kik2){kiq){k2q) 



etiki) 



V2 



Tr 



(1 ± -i^)-fpkiqk2 



(3.9) 



\k,k2){k,q){k2q) 

It is easy to check that this vector satisfies standard conditions for polarization vector: 



0, e^e 



(3.10) 



Then, using formulae of Fierz transformation 



(1 T \{l ± f)7pl,, = 2[1 T ± l%i , 



(l±7Vl,, f(l±7')7, 



(1±7V],J(1±7')7. 



(z, j, fc, / are indices that label the components of 4 x 4-matrices), we obtain: 

(1 T Tr [(1 ± ^'')lpkqk2] = 2(1 T l'')kqk2{l ± 7') = 4(1 T 7')^i#2 

(1 ± 7^)7^ Tr [(1 ± 7')7pA:i^A;2] = -(1 ± 7')7''^igA:2(l ± 7')7p = 
= -2(1 ± j')Ykiqk2lp = 4(1 ± -f')k2qki • 
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Therefore 



and we have 



Y Tr (1 ± 7^)7p^i#2 = 2(1 T l^)kqk2 + 2(1 ± 7')^2gfci 



4J(/ci/c2)(fcig)(M) 



(3.11) 



Note that k2 is an arbitrary massless vector of problem. It follows from this that for 
photon with four-momentum ^2 polarization vectors may be chosen in form ()3.9p too (see 
also jSj): 



Tr 



1 ± -i^)^pkiqk2 



Tr 



'I ip 'y^)jpk2qki 



{kik2){kiq){k2q) 



{kik2){kiq){k2q) 



e^ik2 



(3.12) 



This circumstance simplifies calculations for two-photon process essentially. For example, 
take place expression 



etik^)ef{k2) 



Tr 



(1 ± 7^)7^^171.^2 



Kkik2 



(3.13) 



et al. 

In conclusion of this section we consider one more problem. In reality, polarization 
vectors for the photon with four-momentum ki, satisfied ()3.in|) . may be constructed with 
help of any vectors a and b: 



Tr 



1 ± 7^)7^^106 



ki a b 
ki a b 



4y2G 

Using identity ()2.6|) we have 

Tr[(l ± ^^hpkiab] Tr[(l t -i^)kidckicd\ = Tr[(l ± -f^)-fpkiab]8G 
= Tr[(l T 7^)kiab-fp{l T l^)kidchcd] = Tr[fcia6(l ± 7^)7^^:1^0(1 ± j^)kicd] 
= 4A;ipTr[(l t i^)hdbdckicd] - Tr[fcia6(l ± -i^)kilpdc{l ± 7^)fciCfi] = 
= 4A;ipTr[(l t j^)kidbdckicd] - Tr[(l ± 7^)7pA;icd] Tr[(l ± -f^)kicdkidb] . 



(3.14) 



ki c d 
ki c d 
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Thus 

ih; a, h) = -e, {h; c, ^) , 

y \ki a J \ki c a 

Tt\(1 ^ j^)kiabdckicd\ 
+kip '-^ . 



ki a b J \ki c d 

Since last term may be neglected because of gauge invariance, we obtain that replacement 
of the vectors a, b in polarization vector by another vectors c, d leads to the phase factor 

Tt[{1 ± -f^)kicdkiab] / ± 

, = = — [e [ki] a, b)e^[ki; c, d) 

^ l^fki a b\^fki c d\ 

4 Some remarks about other methods of amplitude 
calculation 

Formula (jl.lOp is a particular case of general equation: 

_ UiZuf TT{QuiUiZufUf) TT{QuiUiZufUf) 

UfQui ~ [ufQui) ■ \ = — — — = ^ _ ' = (4.1) 

\uiZuf\ \uiZuf\ JTi^ZuiU^ZufUf) 



where Z is matrix operator, Z = 7°Z'^7°. Formula 1)4.11) describes all methods of am- 
plitude calculation where Z may be 1 , 7^ , 7° , 1 + 7° et al. Detailed classification of 
different methods is given in j2|- Many methods are analyzed in the book jU] also. 

The most of the methods give more simple expressions for amplitudes in comparison 
with ()1.10|) . For example, if Z = 1 (seef?]), we have: 

-^g^. ^ Tr [{pf + m/)(l + -f%)Qipi + m,)(l + j^hj)] 
A^[m,mf + {piPf)] [1 - {n,nf)] + (p,n/)(p/n,) 

However, the scheme considered here has two essential disadvantages in general case: 

1. There is denominator in ()4.ip . hence the ambiguity of the type - can appear during 
the calculations. For example, in ()4.2)1 the denominator is equal to zero for 

nj - -Ui H — -{Pi H pf) . 

mirrif + [PiPf) ruf 
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2. Each of the fermion hnes obtain phase factor [see (|4.H) ]: 

UiZuf 
\uiZuf\ 

If we have the interference diagrams, this circumstance leads to the fact that expres- 
sions for the amphtudes corresponding to different channels obtain different phase 
factors in general case. In this situation the additional calculations are necessary: 
calculation of the relative phase for interference diagrams, the Fierz transformations 
et al. (see |2]). 

At the same time, the presented method, as mentioned above, has no these disadvantages. 
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